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Abstract 

The  radical  pair  model  of  the  avian  magnetoreceptor  relies  on  long-lived  electron  spin  coherence. 
Dephasing,  resulting  from  interactions  of  the  spins  with  their  fluctuating  environment,  is  generally 
assumed  to  degrade  the  sensitivity  of  this  compass  to  the  direction  of  the  Earth's  magnetic  field.  Here 
we  argue  that  certain  spin  relaxation  mechanisms  can  enhance  its  performance.  We  focus  on  the 
flavin-tryptophan  radical  pair  in  cryptochrome,  currently  the  only  candidate  magnetoreceptor 
molecule.  Correlation  functions  for  fluctuations  in  the  distance  between  the  two  radicals  in 
Arabidopsis  thaliana  cryptochrome  1  were  obtained  from  molecular  dynamics  (MD)  simulations  and 
used  to  calculate  the  spin  relaxation  caused  by  modulation  of  the  exchange  and  dipolar  interactions. 
We  find  that  intermediate  spin  relaxation  rates  afford  substantial  enhancements  in  the  sensitivity  of 
the  reaction  yields  to  an  Earth-strength  magnetic  field.  Supported  by  calculations  using  toy  radical  pair 
models,  we  argue  that  these  enhancements  could  be  consistent  with  the  molecular  dynamics  and 
magnetic  interactions  in  avian  cryptochromes. 


1.  Introduction 

The  rates  and  yields  of  free  radical  reactions  are  often  sensitive  to  applied  magnetic  fields,  a  phenomenon  known 
as  the  radical  pair  mechanism  [1,2].  It  has  been  suggested  [3]  that  migratory  birds  sense  the  direction  of  the 
Earth's  magnetic  field  [4—7]  by  means  of  radical  pairs  formed  photochemically  in  cryptochrome  proteins  located 
in  their  retinas  [8].  Experimental  and  theoretical  support  for  a  radical  pair  mechanism  of  compass 
magnetoreception  is  accumulating  (reviewed  in  [9-17])  although  it  is  not  yet  clear  whether  cryptochrome  is 
actually  the  sensor. 

In  vitro  photo -excitation  of  the  flavin  adenine  dinucleotide  (FAD)  cofactor  in  cryptochromes  leads  to  the 
formation  of  radical  pairs  via  sequential  electron  transfers  along  the  ‘Trp-triad’,  a  chain  of  three  tryptophan 
residues  within  the  protein  [18-21].  Donation  of  an  electron  to  the  excited  singlet  state  of  the  FAD  leads  to  the 
flavin  radical,  FAD*-,  and  the  radical  form  of  the  terminal  tryptophan,  TrpHc’+  (supplementary  figure  SI).  This 
[FAD*-  TrpHc ,+]  radical  pair  in  both  Arabidopsis  thaliana  cryptochrome- 1  (AtCryl)  and  a  closely  related  DNA 
photolyase  has  been  shown  to  respond  to  weak  magnetic  fields  (~  1  mT)  and  has  at  least  some  of  the  properties 
required  of  a  magnetic  compass  sensor  [22, 23].  The  origin  of  the  magnetic  field  effect  lies  in  a  combination  of 
the  spin  dynamics  of  the  radical  pair  and  its  spin- selective  recombination  reactions,  figures  1  (a)  and  (b)  and 
supplementary  figure  S 1  [9] . 

Various  spin  interactions  determine  the  magnetic  field  effect.  The  hyperfine  interactions  [24]  of  the  electron 
spins  with  magnetic  nuclei  (mainly  and  14N)  result  in  coherent,  oscillatory  interconversion  of  the  electronic 
singlet  and  triplet  states  of  the  radical  pair.  The  electronic  Zeeman  interaction  with  the  external  magnetic  field 
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(b) 


(c) 

S  ^  T  RP0 


Figure  1 .  Radical  pair  reaction  schemes.  S  and  T  are  the  singlet  and  triplet  states  of  the  radical  pair.  The  curved  orange  arrows  indicate 
the  coherent  interconversion  of  the  S  and  T  states  by  hyperfine  and  Zeeman  interactions.  ks ,  kT ,  kh,  and  kf  are  rate  constants.  The 
electron  Zeeman  interaction  alters  the  extent  and  timing  of  S  T  interconversion  in  the  radical  pair  so  that  the  final  yields  of  the 

reaction  products,  once  all  the  radicals  have  reacted,  depend  on  the  orientation  of  the  radical  pair  with  respect  to  the  magnetic  field,  (a) 
A  generic  radical  pair  reaction  scheme  in  which  the  S  and  T  states  are  able  to  react  spin-selectively  to  form  singlet  and  triplet  products, 
Ps  and  PT,  respectively,  (b)  A  reaction  scheme  for  [FAD’-  TrpHc'+]  (RPi)  in  AtCryl  in  which  triplet  recombination  cannot  occur. 
Instead,  the  S  and  T  states  undergo  a  non-spin-selective  reaction  (protonation  of  FAD’-,  deprotonation  of  TrpHc‘+,  or  both)  with  the 
same  rate  constant  to  form  a  product  PF.  (c)  An  extension  of  (b)  in  which  RP0  ([FAD’-  TrpHB'+  ])  is  formed  reversibly  by  electron 
hopping  (TrpHB  TrpHc'+  ^2  TrpHB'+  TrpHc)  during  the  lifetime  of  RPi.  In  both  (b)  and  (c),  the  singlet  product,  Ps,  is  the  ground 
state  of  the  protein. 


modifies  this  time- dependence,  altering  the  probability  of  spin- dependent  radical  recombination  and  changing 
the  yields  of  the  reaction  products  (figures  1  (a)  and  (b)).  The  anisotropy  of  the  hyperfine  interactions  causes  the 
effect  of  the  external  field  to  depend  on  its  direction,  as  required  for  a  compass  sensor.  Two  other  important 
interactions — usually  neglected  in  this  context — are  the  electron-electron  exchange  and  dipolar  couplings 
[25, 26],  both  of  which  depend  on  r,  the  distance  between  the  centres  of  electron  spin  density  in  the  two  radicals. 
The  exchange  interaction  is  usually  assumed  to  have  the  exponential  form: 

J(r)  =  /0exp(-/3r),  (1) 

with  a  range  parameter  /3  =  14  n  m  1  for  radicals  embedded  in  a  protein  [25, 27].  The  strength  of  the 
interaction,  /0,  has  a  complex  dependence  on  the  electronic  properties  of  the  radicals  and  is  difficult  to  calculate 
with  any  accuracy.  The  dipolar  interaction,  however,  can  reliably  be  obtained  by  assuming  the  electron  spins  to 
be  point  dipoles  (normally  an  excellent  approximation  for  radical  pairs  [  ]): 

D(r)  =  -2-785mT  =  -77-80MHz.  (2) 

(r/nm)3  (r/nm)3 

In  the  absence  of  other  spin  interactions,  /(r)  and  D(r )  remove  the  degeneracy  of  the  singlet  (S)  and  the  three 
triplet  (T+1,  T0,  T_:)  states: 

E( S)  =  +/(r),  E (Tm)  =  -J(r)  +  D(r)(m2  -  0  (3) 

(m,  the  magnetic  quantum  number  of  the  triplet,  equals  0  or  ±  1).  These  interactions  hinder  efficient  singlet- 
triplet  interconversion  and  so  degrade  the  sensitivity  of  the  radical  pair  to  the  geomagnetic  field  [25]. 

The  magnetically  sensitive  [FAD*-  TrpHc’+]  radical  pair  in  AtCryl  [23]hasr  =  1.90nm[28], 
giving  D(r)  =  — 11  MHz  and  1/  (r)  I  =  6.3  MHz  (assuming  l/0 1  «  8.0  x  107Tand  f3  =  14nm-1[25]). 
Experimental  measurements  of /(r)  for  [FAD*-  TrpHc’+  ]  in  two  DASH-type  cryptochromes  are  in  rough 
agreement  with  this  estimate:  +2.0  MHz  for  Xenopus  laevis  [29]  and  +4.2  MHz  for  Synechocystis  [30].  A  more 
exact  calculation  of  D(r),  avoiding  the  point-dipole  approximation,  confirmed  the  —  1 1  MHz  value  to  within 
0.6  MHz  [25].  It  is  therefore  clear  that  both  1/  (r)  I  and  ID  (r)  I  are  likely  to  be  larger  than  the  electron  Zeeman 
interaction  (~1.4  MHz)  with  the  Earth's  magnetic  field  (~50  //T),  and  cannot  therefore  be  dismissed. 

Although  values  of  1/  (r)  I  or  ID  (r)  I  larger  than  a  few  MHz  would  be  sufficient  to  mar  the  performance  of  a 
radical  pair  as  a  geomagnetic  sensor,  the  effects  of  the  two  interactions  can  partially  cancel  under  certain 
circumstances,  specifically  when  D(r)//(r)  «  +6  or  —3  [25].  Using  estimates  of  l/0 1  derived  from  data  on 
photosynthetic  reaction  centre  proteins  [3 1-33],  it  has  been  argued  that  such  ratios  can  be  expected  when 
r  «  2.0  ±  0.2  nm,  i.e.  a  separation  close  to  that  of  FAD*-  and  TrpHc*+  in  cryptochromes  (~1.9  nm)  [2  ]. 
However,  the  values  of  /(r)  and  D(r)  are  not  known  for  any  of  the  four  avian  cryptochromes,  and  it  has  not  been 
established  how  effective  the  cancellation  would  be  for  radicals  as  complex  as  flavin  and  tryptophan.  It  is 
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therefore  appropriate  to  ask  whether  the  damaging  effects  of  /(r)  and  D(r)  could  be  circumvented  in  some  other 
way.  One  possibility,  discussed  later,  is  that  £the  quantum  Zeno  effect  immunises  the  avian  compass  against  the 
deleterious  effects  of  exchange  and  dipolar  interactions’  [34].  Here,  we  examine  whether  electron  spin  relaxation 
could  come  to  the  rescue. 

Spin  relaxation,  the  process  that  drives  the  electron-nuclear  spin  system  of  the  radical  pair  towards  thermal 
equilibrium  [24],  is  normally  expected  to  attenuate  the  sensitivity  of  the  compass  by  destroying  the  spin 
coherence  that  is  essential  for  its  operation  [35].  It  seems  plausible  that  the  optimum  performance  would  require 
negligible  loss  of  coherence  during  the  lifetime  of  the  radical  pair.  There  have  been  several  theoretical  treatments 
of  spin  relaxation  in  the  context  of  magnetoreception,  almost  all  of  which  have  dealt  with  generic  relaxation 
mechanisms  furnished  with  phenomenological  rate  constants  [36-41].  Intriguingly,  a  few  of  these  studies  have 
indicated  that  the  compass  sensitivity  could  be  enhanced  by  the  inclusion  of  certain  abstract  spin  relaxation 
pathways  [37, 38, 40, 42, 43].  In  the  following,  we  discuss  the  spin  relaxation  that  arises  from  realistic  spin 
interactions  modulated  by  realistic  protein  dynamics.  Specifically,  we  show  that  the  time  dependence  of  /(r)  and 
D(r)  arising  from  stochastic  fluctuations  in  the  separation  of  the  FAD*-  and  TrpHc*+  radicals  in  cryptochrome 
has  the  potential  to  boost  the  magnetic  field  effect  and  so  counteract  the  detrimental  effects  these  interactions 
would  have  in  the  absence  of  molecular  motion. 


2.  /-modulation  by  fluctuations  in  r 


We  start  by  establishing  the  principle  that  the  harmful  effect  of  an  exchange  interaction  on  a  radical  pair  compass 
can  be  alleviated  by  fluctuations  in  the  strength  of  that  interaction.  To  do  so,  we  chose  the  simplest  possible 
radical  pair  comprising  two  electron  spins,  one  of  which  is  coupled  to  a  single  spin-  ^  nucleus  via  an  anisotropic 
hyperfine  interaction.  The  eigenvalues  of  the  anisotropic  hyperfine  tensor  [24],  A^  =  Ayy  =  0, 

Azz  =  42  MHz,  were  chosen  to  resemble  the  hyperfine  tensor  of  the  N5  interaction  in  FAD*-  [44].  Singlet  and 
triplet  pairs  were  allowed  to  react  to  form  distinct  products  (figure  1  (a))  with  identical  rate  constants, 
k  =  ks  =  kT  =  106  s-1,  such  that  the  radical  pair  lifetime  was  1  /is.  At  time  t  =  0,  the  radical  pair  was  in  a  pure 
singlet  state  with  no  nuclear  polarisation.  The  performance  of  this  toy  compass  was  judged  by  the  anisotropy  of 
the  singlet  yield: 

AT>S  =  T>s(0)  -  T>s(90°),  (4) 

where  T>s  (6)  is  the  final  yield  of  the  singlet  product  (Ps  in  figure  1  (a))  and  0  is  the  angle  between  the  external 
magnetic  field  (50  /IT)  and  the  symmetry  axis  of  the  hyperfine  tensor. 

<FS  (0)  was  determined  using  the  equation  of  motion  for  the  radical  pair  density  operator,  p  ( t ): 

^42  =  —  i [h0,  pm  -  kp(t )  -  kSTD(Psp(t)PT  +  pTp(t)Ps)  =  -hp(t).  (5) 

at 

H0  is  the  time-independent  part  of  the  spin  Hamiltonian: 

Ho  =  7e^o [ (^Ax  +  5b*) sin  0  +  (SAz  +  5bz)cos  6]  +  Azz Saz Iz  —  /av^2SA  •  5b  +  — (6) 


in  which  the  three  terms  are,  respectively,  the  electron  Zeeman,  electron-nuclear  hyperfine,  and  average 
exchange  interactions.  SA  and  5B  are  the  electron  spin  angular  momentum  operators  of  the  two  radicals,  and  I  is 
the  spin  operator  of  the  nucleus  in  radical  A.  B0  =  50  /IT  is  the  strength  of  the  external  magnetic  field.  The  other 
terms  in  equation  (5)  account  for  the  radical  recombination  reactions  and  the  effects  of  fluctuations  in  /(r).  The 
latter  is  derived  in  appendix  A.  k$  td  is  the  rate  constant  of  the  spin  relaxation  caused  by  /-modulation,  and  P 
and  P  are  the  singlet  and  triplet  projection  operators.  If  the  fluctuations  in  the  exchange  interaction  occur  on  a 
faster  timescale  than  the  coherent  spin  dynamics,  k$ td  is  given  by  (appendix  A) 

*:std  =  4<[/(r(t))  -/av]2>rc,  (7) 


where  ( •  •  • )  indicates  a  time  average,  /av  =  (/  (r  ( t )) ) ,  and  the  correlation  time  rc  characterises  the  timescale  of  the 
fluctuations  in  r.  The  rate  of  relaxation  is  thus  proportional  to  the  mean  square  deviation  of  the  exchange 
parameter  from  its  mean  value  and  is  faster  the  slower  the  fluctuations.  It  is  clear  from  equation  (5)  that  /- 
modulation  relaxes  density  matrix  elements  of  the  type  ( S  |  p{t)  \  Tm) ;  it  is  therefore  referred  to  as  a  'singlet- 
triplet  dephasing'  (STD)  mechanism  [45].  We  note  that  the  relaxation  term  in  equation  (5)  can  also  be  obtained 
phenomenologically  using  Lindblad  noise  operators,  Lj  [36, 46]: 


dp(t) 
d  t 


+  i[fTo>  P(t)]  +  kp(t)  —  ksjD 


t(£/4p(o  +  pmAi) 


(8) 


withlj  =  {PS,  PT}or2  1/2(PS  —  PT ). 
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(a) 


(c) 


(b) 


(d) 


log,o(kSTo/  s'1),  log10(/cR:R  /  s') 


Figure  2.  Anisotropy  of  the  singlet  yield  ( A4>s)  of  a  toy  radical  pair  as  a  function  of  the  spin  relaxation  rate  constant.  The  solid  lines  in 
(a)  and  (b)  show  A4>s  (equation  (4))  for  the  STD  and  RFR  mechanisms,  respectively,  for  four  values  of  Jav.  The  dashed  lines  are  A4>s  in 
the  absence  of  spin  relaxation,  (c)  The  singlet  yield,  4>s  ( 6 ),  when  the  magnetic  field  is  parallel  ( 6  =  0)  and  perpendicular  ( 0  =  90°)  to 
the  symmetry  axis  of  the  hyperfine  interaction,  for  the  STD  and  RFR  mechanisms  with  Jav  =100  p T.  (d)  A4>s  for  the  STD  relaxation 
mechanism  induced  by  electron  hopping.  Blue  dots:  A4>s  (ho)  calculated  with  explicit  two-site  electron  hopping;  orange  line: 

AT>S  (^std)  calculated  using  equation  (9).  The  green  line  indicates  the  value  of  ki0  ( 1 . 1  x  1 07  s“ *)  determined  from  experiments  on 
AtCryl  [23].  (a)-(c)D  =  0,  magnetic  field  =  50  pT,  principal  components  of  the  hyperfine  interaction,  (A**,  Ayy,  Azz)  =  (0.0, 0.0, 

+  1.5)  mT,  recombination  rate  constants  =  ks  =  kT  =  106s_1  (figure  1(a)).  (d)  As  for  (a)-(c),  with /(RP0)  =  800  mT, 

/(RP)i  =  0.1  mT,  fc0i  =  1.0  x  10los_1,  kh  —  4.9  x  105  s_1,  h  =  2.5  x  105  s_1  (see  figure  1(b)  for  definitions  of  the  rate  constants). 


Finally,  the  singlet  yield  once  all  radicals  have  recombined  is: 

J't  OO  q  r*  OO  q  ^ 

Tr[F p(t)]dt=  k  Tr[Psexp(-fit)p(0)]df,  (9) 

o  Jo 

£  1  ~s 

where  the  superoperator  D  is  defined  in  equation  (5)  and  p  (0)  =  -P  . 

Figure  2(a)  plots  the  anisotropy  A<FS  (equation  (4))  versus  fcSpD  for  four  values  of  /av.  When  fcSTD  is  small 
(e.g.  103  s-1),  there  is  essentially  no  relaxation  during  the  1  /is  lifetime  of  the  radicals.  In  this  limit,  A<FS  is 
dramatically  reduced  when  the  average  exchange  interaction  is  comparable  to  or  stronger  than  the  geomagnetic 
field  (50  /iT).  However,  when  /av  ^  0,  values  of  ks td  greater  than  ~106  s-1  cause  A<FS  to  rise  towards  the  values 
appropriate  for  /av  =  0.  When  the  relaxation  is  faster  still  (/cSTD  >  108  s_1),  A^>s  drops  towards  zero. 
Fluctuations  in /(r)  resulting  in  intermediate  relaxation  rates  (106s_1  <  /cSTD  <  108  s_1)  therefore  boost  the 
compass  sensitivity  well  above  the  level  expected  for  a  time- independent  exchange  interaction  (dashed  lines  in 
figure  2(a)).  Similar  behaviour  was  found  for  modulation  of  the  dipolar  interaction  (appendix  B  and 
supplementary  figure  S2)  and  for  the  reaction  scheme  in  figure  1  (b)  (supplementary  figure  S3). 

To  put  figure  2(a)  in  context,  /-modulation  was  compared  with  a  ‘random  fields’  relaxation  (RFR) 
mechanism  in  which  unspecified,  uncorrelated  molecular  motions  cause  the  expectation  values  of  all  three  spin 
angular  momentum  operators  ( Sx ,  Sy,  and  Sz)  of  both  electrons  to  relax  independently  with  the  same  rate 
constant,  kRFR  (appendix  C).  As  shown  in  figure  2(b),  this  relaxation  mechanism  does  not  enhance  A<FS  but 
causes  it  to  fall  towards  zero  at  much  smaller  values  of  the  rate  constant  than  was  the  case  in  figure  2(a).  Clearly, 
not  all  relaxation  mechanisms  are  able  to  rescue  the  sensitivity  of  a  radical  pair  compass. 

To  shed  light  on  the  differences  between  figures  2(a)  and  (b),  figure  2(c)  plots  <FS  ( 6 )  for  9  =  0  and  9  =  90° 
as  a  function  of  the  relaxation  rate  constant  for  the  two  mechanisms  (with  /av  =  100  //T).  It  can  be  seen  that  very 
fast  relaxation  (>10los-1)  destroys  the  anisotropy  either  by  locking  the  radical  pair  into  its  initial  singlet  state 
(STD)  or  by  relaxing  the  spins  to  the  equilibrium  state  in  which  the  singlet  probability  is  0.25  (RFR).  For  rate 
constants  between  106  s-1  and  108  s-1,  STD  is  much  less  effective  at  relaxing  T>s  (6)  towards  0.25  when  the 
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magnetic  field  is  parallel  to  the  hyperfine  axis  (i 9  =  0)  than  when  it  is  perpendicular  (0  =  90°).  For  the  STD 
mechanism,  therefore,  the  anisotropy  in  the  reaction  yield  is  considerably  larger  when 

106  S_1  <  Zcstd  ^5  108  S_1  than  in  the  absence  of  relaxation  (/cSTD  <  104  s_1).  The  RFR  mechanism,  by  contrast, 
relaxes  the  spins  isotropically  and  so  drives  AT>S  monotonically  to  zero  as  /cRFR  is  increased. 

3.  /- modulation  by  electron  hopping 

STD  has  been  invoked  to  account  for  the  observed  magnetic- field  dependence  of  the  flavin-tryptophan  radical 
pair  reaction  in  AtCryl  [23].  /-modulation  was  considered  to  arise  from  reversible  electron  hopping  between  the 
distal  and  medial  members  of  the  Trp-triad,  i.e.  [FAD‘-TrpHB  TrpHc'+]  [FAD*-  TrpHB  +  TrpHc] 

(figure  1  (c)).  As  described  above,  the  former  (denoted  RPO  has  rx  =  1 .90  nm  and  exchange  and  dipolar 
interactions  of  a  few  MHz;  the  latter  (denoted  RP0)  has  a  centre -to -centre  distance  r0  =  1.32  nm  in  AtCryl  and 
estimated  interactions:  I /  (r0)  1^21  GHz,  D(r0)  =  —34  MHz  [25].  With  such  a  large  exchange  interaction  in 
RP0,  reversible  electron  hopping  results  in  strong  modulation  of  /  (r)  and  potentially  rapid  relaxation  of  singlet- 
triplet  coherences  in  RP:.  Provided  fc01,  the  electron  hopping  rate  constant  for  RP0  — ►  RP1?  is  not  much  larger 
than  I /  (r0)  I ,  then  every  time  RP:  is  converted  into  RP0  the  coherence  experiences  a  large  phase  shift  before  the 
electron  jumps  back.  Under  these  conditions,  the  dephasing  rate  constant  ks td  in  equation  (5)  is  expected  to  be 
equal  to  Zc10,  the  rate  constant  for  RP0  <—  RP:  [23]. 

This  prediction  can  be  verified  with  a  toy  model  similar  to  that  used  above.  Figure  2(d)  shows  the  singlet  yield 
anisotropy  determined  (i)  as  a  function  of  ks td  using  equation  (5)  and  (ii)  as  a  function  of  Zc10  by  solving  a 
Liouville  equation  containing  explicit  reversible  electron  hopping  between  RP0  and  RPl5  coherent  spin 
dynamics  in  both  radical  pairs,  and  spin -selective  recombination  for  RP:  (figure  1  (c)  and  appendix  D).  The 
parameters  for  method  (ii)  were  chosen  to  approximate  those  determined  from  the  observed  magnetic  field 
effects  on  [FAD*-  TrpHc'+]  in  AtCryl  [23].  The  two  calculations  are  in  excellent  agreement  confirming  that  the 
dephasing  rate  constant  is  indeed  given  by  Zc10.  Note  also  that  the  value  of  fc10  (=  1-1  x  107  s-1)  determined  from 
the  magnetic  field  effects  on  AtCryl  [23]  corresponds  to  a  position  in  figure  2(d)  close  to  the  maximum 
enhancement  of  A<f>s .  However,  this  should  be  treated  with  caution  because  of  our  use  of  a  grossly  simplified 
spin  system  for  figure  2(d)  (one  nuclear  spin  instead  of  ~10;  eight  spin  states  instead  of  ~212). 

4.  /- modulation  in  cryptochrome 

We  now  turn  to  STD  spin  relaxation  for  the  [FAD*-  TrpHc ,+]  radical  pair  in  cryptochrome.  Fluctuations  in  the 
separation  and  relative  orientation  of  the  two  radicals  were  characterised  by  means  of  all- atom  molecular 
dynamics  (MD)  simulations  [47].  The  calculations  were  performed  using  a  modified  version  of  the  crystal 
structure  of  AtCryl  (PDB  ID  1U3C  [28])  containing  the  flavin  radical,  FAD*-,  and  the  radical  form  of  the 
terminal  residue  of  the  Trp-triad,  TrpHc*+  (W324)  with  force-field  parameters  and  atomic  charges  as  in  [48]. 
The  protein  was  solvated  in  an  aqueous  solution  of  50  mM  NaCl  in  a  periodic  box  of  dimensions 
103  x  111  x  99  A3  containing  a  total  of  113  455  atoms.  After  extensive  equilibration,  seven  statistically 
independent  MD  trajectories,  three  spanning  300  ns  and  four  covering  100  ns,  were  accumulated  for  the 
microcanonical  ensemble  at  300  K.  The  initial  configuration  of  the  system  was  identical  in  the  seven  simulations; 
stochastic  behaviour  emerged  from  the  Langevin  thermostat  [49],  as  implemented  in  NAMD  [50].  The  protein 
dynamics  were  sampled  at  intervals  of  0.5  or  5  ps.  The  accumulated  trajectories  cover  1.3  /is.  Because  the 
motions  of  the  radicals  are  considerably  faster  than  both  spin  relaxation  and  radical  recombination,  it  is  not 
necessary  for  the  simulations  to  include  the  reactions  or  to  span  the  lifetime  of  the  radicals.  Electron  hopping 
between  TrpHB  and  TrpHc  [23]  was  not  included.  Further  details  of  the  simulations  can  be  found  in  [47]. 

Motional  correlation  functions,  calculated  from  the  MD  trajectories,  were  used  to  obtain  spectral  densities 
and  hence  a  Redfield  relaxation  superoperator  [51, 52]  which,  when  combined  with  (i)  electron  Zeeman  and 
hyperfine  interactions,  (ii)  average  exchange  and  dipolar  interactions,  and  (iii)  spin- selective  recombination 
reactions,  allowed  calculation  of  the  singlet  yield  <4>s  (0>  </>),  where  the  polar  angles  6  and  (j)  specify  the  direction  of 
the  external  50  /llT  magnetic  field  (appendix  A,  equation  (20)).  The  slow  relaxation  arising  from  cross¬ 
correlation  of  the  exchange  and  the  dipolar  interactions  was  neglected.  The  reaction  yield  anisotropy,  A<FS,  was 
calculated  as  the  difference  between  the  maximum  and  minimum  values  of  <4>s  (0>  0).  Spin  relaxation  induced 
by  fluctuations  in  the  hyperfine  interactions  was  the  subject  of  [47]  and  is  not  treated  here.  The  forward  electron 
transfer  steps  along  the  Trp-triad  are  complete  within  a  nanosecond  which  is  too  fast  for  singlet-triplet 
coherence  to  be  generated  in  the  intermediate  radical  pairs  [FAD*-  TrpHA ,+]  and  [FAD*-  TrpHB'+  ]. 
Consequently,  STD  relaxation  is  expected  to  occur  only  in  [FAD*-  TrpHc'+]. 

Figures  3(a)-(c)  show  the  variation  in  r,  /(r)  and  D(r)  over  the  course  of  MD  trajectories  spanning  1 .3  //s  in 
total. /(r)  and  D(r)  were  calculated  using  equations  (1)  and  (2),  with  l/0 1  =  8  x  107Tand  /3  =  14nm-1.The 
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Figure  3.  Variation  in  (a)  the  radical  separation,  (b)  the  exchange  interaction  and  (c)  the  dipolar  interaction  of  [FAD'-  TrpHc'+]  in 
AtCryl .  The  MD  trajectories  were  sampled  every  5  ps.  The  different  colours  indicate  different  MD  runs  starting  from  the  same 
configuration.  The  horizontal  lines  in  the  three  panels  give  the  mean  values:  rav,  /av,  and  Dav.  Shown  on  the  right-hand  side  of  (a)-(c) 
are  the  corresponding  distribution  functions,  (d)  Auto -correlation  functions,  gK  (r)  (equation  (10)),  for  the  exchange  (pink)  and 
dipolar  (green)  interactions  in  [FAD*-  TrpH*+]  in  AtCryl ,  calculated  from  the  MD  trajectories.  The  dashed  black  lines  are  the  multi¬ 
exponential  fits  used  to  calculate  the  relaxation  superoperator;  both  have  been  extrapolated  back  to  r  =  100  fs.  gr  (r)  decays  to  ~0.79 
within  the  first  500  fs;  motions  on  this  timescale  are  far  too  fast  to  cause  spin  relaxation. 


Figure  4.  Molecular  dynamics  structures  of  AtCryl  in  the  region  of  the  FAD  and  Trp  triad.  The  positions  of  the  two  components  of  the 
radical  pair  (FAD’-  and  TrpHc*+ ),  the  proximal  and  medial  tryptophans  (TrpHA  and  TrpHB),  and  the  backbone  in  the  average 
structure  (rav  =  1 .89  nm)  are  shown  in  green  and  grey.  The  positions  of  the  two  radicals,  and  the  backbone  in  the  region  of  TrpHc*+, 
are  shown  in  red  and  blue  for  the  structures  with,  respectively,  the  largest  (rmax  =  2.43  nm)  and  smallest  (rmin  =  1.69  nm)  radical 
separations.  The  curved  arrow  indicates  the  range  of  motion  of  TrpHc'+. 


mean  distance  between  the  centres  of  spin  density  (=b  standard  deviation;  see  supplementary  figure  S4)  is 
rav  =  1.89  zb  0.05  nm.  The  mean  interaction  strengths  (=b  standard  deviation)  are  l/av  I  =  8.42  db  5.64  MHz 
andDav  =  —11.5  =b  0.9  MHz.  The  maximum  and  minimum  separations  are  rmax  =  2.43  nm  and 
rmin  =  1.69  nm,  so  that  l/(rmin)l  =  115  MHz,  l/(rmax)l  =  3.85  kHz,D(rmin)  =  —16.1  MHz,  and 
D(rm  ax)  =  -5.5  MHz. 

Figure  4  gives  an  impression  of  the  positions  of  the  two  radicals  in  AfCry  1  in  the  average  structure  and  in  the 
two  extreme  conformations,  r  =  rmin  and  r  =  rmax.  The  variation  in  the  radical  separation  is  mainly  attributable 
to  the  greater  mobility  of  the  TrpHc '+  radical  compared  to  the  more  snugly  bound  FAD’-  radical  located  in  the 
centre  of  the  protein.  Representations  of  the  semi-occupied  molecular  orbitals  and  the  spin  densities  in  the  two 
radicals  are  shown  in  supplementary  figure  S4. 
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Figure  5.  Dependence  of  (a)  the  average  singlet  yield  (<f>s)  and  (b)  the  fractional  anisotropy  of  the  singlet  yield  (r  =  AT>s  /  4>s) on  the 
mean  exchange  interaction,  /av.  (c)  Dependence  of  r  on  the  radical  pair  lifetime  for  selected  values  of  /av.  The  green  and  orange 
lines  are  for  radical  pairs  respectively  with  and  without  the  spin  relaxation  caused  by  fluctuations  in  the  exchange  and  dipolar 
interactions.  For  all  three  panels,  the  magnetic  field  was  50  /iT.  In  (a)  and  (b)  the  radical  pair  lifetime  was  r  =  fcs-1  =  h:~l  =  1  /^s. 
[FAD’-  TrpHc'+]  was  modelled  by  including  the  N5  and  N10  hyperfine  interactions  of  FAD’-  in  one  radical  and  no  hyperfine 
interactions  in  the  other  radical  (i.e.  a  total  of  36  spin  states).  The  hyperfine  and  average  electron-electron  dipolar  interaction  tensors 
are  given  in  the  supporting  information  of  [47].  The  relaxation  effects  in  the  lower  two  panels  in  (c)  arise  entirely  (left)  and  mostly 
(right)  from  D-modulation. 


Given  the  uncertainty  in  the  value  of  /0  (equation  (1)),  we  henceforth  treat  I /av  I  as  an  adjustable  parameter 
and  model  the  exchange  interaction  as  /(r)  =  /avexp(— (3(r  —  rav)),  with/5  =  14  nm”1  andrav  =  1.89  nm  as 
above.  Except  for  very  small  values  of  l/av  I ,  the  fluctuations  in  /  (r)  are  much  larger  than  those  in  D  (r),  meaning 
that  /-modulation  should  be  the  dominant  source  of  spin  relaxation. 

The  auto-correlation  functions 

gK(r)  =  (l K(t )  -  Kw][K(t  +  r)  -  Kw])  (10) 

( K  =  r,  /,  or  D,  and  Kav  =  (K  ( t ) ) )  that  characterise  the  relevant  motions  of  the  radical  pair  are  highly  non¬ 
exponential  as  a  result  of  motions  on  timescales  ranging  from  picoseconds  to  nanoseconds  (figure  3(d)).  The 
average  time  constant  for  the  decay  of  gr  (r)  (8.5  ns)  is  determined  mainly  by  low- amplitude  components  in  the 
ranger  ^  10  ns.  In  combination  with  the  large  fluctuations  in  /(r(t))  noted  above,  these  comparatively  slow 
motions  tend  to  dominate  the  spin  relaxation.  By  contrast,  the  rapidly  decaying  components  of  gr  (r ) 

(r  <  500  fs)  are  of  minor  importance.  The  average  time  constants  for  gj  (r)  and  gD  (r )  are  8.5  ns  and  7.6  ns, 
respectively,  reflecting  the  different  r-dependence  of  the  two  interactions. 

Figures  5(a)  and  (b)  show,  respectively,  the  singlet  yield  averaged  over  all  directions  of  the  magnetic  field,  T>s, 
and  the  relative  anisotropy,  T  =  A<Ts  /  5>s>  both  as  functions  of  /av.  The  lifetime  of  the  radical  pair 
(r  =  /cs-1  =  kj 1)  was  1  //s.  When  l/av  I  <  5  MHz,  the  protein  dynamics  have  a  minimal  effect  on  <FS 
(figure  5(a)).  For  larger  l/av  I,  relaxation  tends  to  reduce  4>s  towards  the  statistical  value  of  0.25.  The  effect  of  spin 
relaxation  on  the  anisotropy  (figure  5(b))  is  small  when  I /av  I  is  outside  the  range  2  MHz  <  l/av  I  <  25  MHz  but 
shows  a  complex  behaviour  within  this  range.  For  most  of  these  intermediate  values  of  l/av  I ,  STD  boosts  the 
anisotropy  above  that  found  when  relaxation  is  excluded.  T  is  particularly  strongly  enhanced  for  positive  values 
of  the  exchange  interaction  (7  MHz  <  /av  <  18  MHz). 

Figure  5(c)  offers  additional  insight  by  plotting  the  relative  singlet  yield  anisotropy  T  against  the  radical  pair 
lifetime  for  a  few  representative  values  of  /av.  In  the  absence  of  relaxation  (orange  lines),  T  is  small  for  very  short 
lifetimes  (r  <  10  ns),  grows  as  r  approaches  the  reciprocal  of  the  Farmor  period  (~700  ns),  and  eventually 
reaches  a  plateau  when  r  >  ~  1 0  /is.  When  relaxation  is  included  (green  lines)  and  I /av  I  is  small  (e.g.  0.0  and 
—0.75  MHz),  the  anisotropy  is  always  less  than  or  equal  to  that  in  the  absence  of  relaxation.  For  larger  /av 
however  (e.g.  +16.8  MHz  and  +9.0  MHz),  the  relaxation  causes  a  substantial  increase  in  T  for  radical  pairs  with 
intermediate  lifetimes  (10  ns  <  r  <  1  /is).  For  example,  when /av  =  +16.8  MHz,  and  r  =  100  ns,  the  STD 
mechanism  boosts  T  by  a  factor  of  34.  No  enhancement  of  any  magnitude  was  found  in  our  earlier  study  of  the 
relaxation  resulting  from  fluctuations  in  the  hyperfine  interactions  in  FAD*-  and  TrpHc'+  [47]. 
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In  summary,  unless  the  mean  exchange  interaction  is  too  large  or  too  small,  the  thermal  fluctuations  in  the 
separation  of  the  FAD*-  and  TrpHc'+  radicals  in  cryptochrome  can  induce  spin  relaxation  of  such  a  rate  as  to 
enhance  the  magnetic  sensitivity,  sometimes  quite  substantially. 


5.  Discussion 

In  the  following  paragraphs  we  discuss  other  ways  in  which  Nature  might  have  avoided  the  problems  posed  by 
exchange  and  dipolar  interactions  in  a  [FAD*-  TrpHc'+]  cryptochrome  magnetosensor.  We  then  comment  on 
possible  reasons  why  the  sensitivity  enhancements  afforded  by  fluctuations  in  the  radical  separation  might  not 
be  manifested  in  cryptochrome. 


5.1.  Asymmetric  recombination  rates 

Dellis  and  Kominis  have  suggested  that  harmful  effects  of  exchange  and  dipolar  interactions  maybe  avoided  if 
the  radical  recombination  kinetics  are  highly  asymmetric  [34].  They  argued  that  the  anisotropy  AT>s  can  be 
larger  when  the  rate  constants  for  recombination  of  the  singlet  and  triplet  states  (figure  1  (a))  are  such  that 
kT  ks  rather  than  kT  =  ks.  We  believe  this  is  unlikely  to  occur  in  a  cryptochrome  magnetoreceptor,  for  two 
reasons.  First,  our  experience  with  toy  models  (supplementary  figure  S5)  suggests  that  any  benefit  arising  from 
the  condition  /cT  ks  only  accrues  when  the  average  exchange  and/  or  dipolar  interactions  are  stronger  than 
the  effective  hyperfine  interactions  in  both  radicals,  viz.  39  MHz  =  1.40  mT  for  FAD*-  and  54  MHz  =  1.94  mT 
for  TrpHc  *+  [23] .  A  value  of  l/av  I  in  excess  of 40-50  MHz  would  result  in  such  a  small  singlet  yield  anisotropy 
(see,  e.g.,  figure  5(b))  that  it  would  probably  be  beyond  rescue  by  even  the  most  favourable  kinetics  [53].  Second, 
a  reaction  scheme  in  which  singlet  and  triplet  radical  pair  states  both  recombine  spin  selectively  to  give  distinct 
products  (figure  1(a))  is  not  appropriate  for  [FAD*-  TrpHc*+]  because  there  is  no  energetically  accessible 
molecular  triplet  state  in  cryptochrome  to  which  the  triplet  radical  pair  could  recombine  [23].  Rather,  one 
should  consider  a  scheme  (figure  1  (b))  in  which  the  singlet  recombines  spin  selectively  (to  form  the  ground  state 
of  the  protein,  with  rate  constant  /cb)  in  competition  with  a  reaction  such  as  (de)protonation  of  one  or  both  of  the 
radicals  [23]  whose  rate  constant  (kf)  is  independent  of  the  spin  state  of  the  pair.  For  such  a  reaction  scheme, 
there  is  no  advantage  to  be  gained  from  highly  asymmetric  rate  constants  (supplementary  figure  S5). 


5.2.  Alternative  radical  pairs 

Despite  the  finding  that  [FAD*-  TrpHc ,+]  in  isolated  AtCryl  in  vitro  responds  to  relatively  weak  magnetic  fields 
[23],  it  is  not  clear  that  this  radical  pair  would  necessarily  be  the  sensor  in  an  in  vivo  cryptochrome 
magnetoreceptor.  Protein-protein  and  protein-ligand  interactions  in  a  cellular  environment  could  result  in  a 
different,  and  possibly  more  sensitive,  radical  pair  via  an  alternative  or  extended  electron  transfer  pathway 
[44, 54, 55].  It  is  therefore  important  to  ask  whether  Nature  could  have  contrived  a  cryptochrome-based 
magnetic  sensor  in  which  the  exchange  and  dipolar  interactions  are  so  weak  that  they  do  not  adversely  affect  its 
performance.  This  is  a  complex  question  with  no  clear  answer.  Using  equation  (1)  with  =  14  nm-1,  an 

increase  in  r  of  only  0.33  nm  would  result  in  a  100-fold  reduction  in  /(r),  making  it  negligibly  small.  This  might 
be  achieved  if  a  fourth  tryptophan  residue  were  involved  in  the  electron  transfer  chain  that  creates  the 
magnetically  sensitive  radical  pair  [56, 57] .  However  such  an  increase  in  the  separation  of  the  radicals  would  also 
slow  down  the  reverse  electron  transfer  reaction  by  about  two  orders  of  magnitude,  prolong  the  lifetime  of  the 
radical  pair  and  make  it  much  more  vulnerable  to  loss  of  coherence  by  spin  relaxation.  This  problem  could  be 
avoided  by  multi-step  reverse  electron  transfer  along  the  Trp-triad,  e.g.  [FAD*-  TrpHB  TrpHc'+]  — ►  [FAD*- 
TrpHB+TrpHc]  — ►  FAD  +  TrpHB  +  TrpHc,  although  it  is  not  clear  that  this  would  be  energetically  feasible 
[58, 59].  Such  electron  hopping  could  be  accompanied  by  strong /-modulation  as  discussed  in  section  3. 

Turning  now  to  the  dipolar  interaction,  ID  (r)  I  only  becomes  comparable  to  the  Zeeman  interaction  with  the 
geomagnetic  field  (50  /jlT  =  1.4  MHz)  when  r  =  3.8  nm  (equation  (2)).  In  the  x-ray  structure  of  AtCryl  [28], 
the  tryptophan  and  tyrosine  residues  that  are  furthest  from  the  FAD  are  W447  at  2.90  nm  and  Y1 70  at  2.44  nm 
(edge-to-edge  distances  between  aromatic  rings).  Neither  is  far  enough  away  to  have  ID  (r)  I  <  1.4  MHz  and 
neither  is  part  of  an  obvious  electron  transport  chain  linking  it  to  the  FAD.  One  way  to  realise  a  radical- radical 
separation  exceeding  3.8  nm  would  be  if  the  ultimate  electron  donor  were  located  on  another  (unknown) 
protein  complexed  to  the  cryptochrome.  Once  again,  a  sequential  reverse  electron  transfer  pathway  would  be 
needed  to  allow  fast  enough  singlet  recombination.  Finally,  an  effective  radical  separation  as  large  as  3-4  nm 
could  conceivably  be  contrived  if  the  FAD*-  inside  the  cryptochrome  were  paired  with  a  freely  diffusing  radical 
formed,  for  example,  by  electron  transfer  to  the  Trp  radical  in  [FAD*-  TrpHc ,+  ].  Ascorbic  acid  has  been 
discussed  in  this  context  [44]. 
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In  summary,  although  one  can  imagine  cryptochrome -derived  radical  pair  magnetic  sensors  in  which  /av 
and  Dav  are  small  enough  that  they  do  not  present  a  problem,  such  arrangements  do  not,  at  present,  seem  very 
likely. 

5.3.  Relaxation-enhanced  compass  sensitivity  in  cryptochrome 

The  results  presented  in  figure  5  suggest  that  the  anisotropic  response  of  [FAD*-  TrpHc'+]  in  AtCry  1  is 
particularly  strongly  enhanced  for  positive  /av  in  the  range  ~7  MHz  to  ~  1 8  MHz.  These  values  are  of  the  same 
sign  but  somewhat  larger  than  the  only  experimental  measurements  of  this  quantity  (in  different 
cryptochromes):  2.0  MHz  [29]  and  4.2  MHz  [30].  However,  it  should  be  borne  in  mind  that  the  exchange 
interaction  depends  on  the  separation  of  the  radicals  and  the  surrounding  amino  acid  residues  and  could  be 
larger  in  an  avian  magnetoreceptor  cryptochrome  than  in  a  cryptochrome  with  a  different  biological  function. 

Although  we  have  focussed  here  on  the  modulation  of  exchange  and  dipolar  interactions,  this  is  not  the  only 
source  of  spin  relaxation.  The  other  obvious  contributors  are  (i)  modulation  of  anisotropic  hyperfme 
interactions  by  librational  motions  of  the  radicals  within  their  binding  sites  and  (ii)  modulation  of  isotropic 
hyperfme  interactions  by  conformational  fluctuations,  in  particular  in  dihedral  angles  [47].  Further  work  will  be 
needed  to  determine  (a)  whether  the  enhancements  discussed  here  survive  in  the  presence  of  these  other  spin- 
relaxation  pathways  and  (b)  the  consequences  of  /-modulation  and  electron  hopping  occurring  in  tandem. 

The  simulations  of  [FAD*-  TrpHc’+]  described  here  were  performed  using  a  drastically  simplified  spin 
system,  comprising  two  electron  spins  and  two  nuclear  spins  whose  hyperfme  tensors  are  those  of  the  two  14N 
nuclei  (N5  and  N10)  in  the  central  ring  of  the  FAD*-  isoalloxazine  ring  system  [44].  These  two  interactions — 
strongly  axially  anisotropic  with  parallel  symmetry  axes — dominate  the  spin  dynamics  of  FAD*-  and  arguably 
make  flavin  radicals  near-optimal  as  components  of  a  magnetic  compass  sensor  [44, 60].  However,  the  influence 
of  other  nuclear  spins  cannot  be  discounted.  In  particular,  a  number  of  the  hyperfme  interactions  in  TrpHc’+ 
are  significantly  anisotropic  and  of  lower  symmetry  [44] .  It  seems  likely  that  the  lower  the  overall  anisotropy  of 
the  radical  pair,  the  less  scope  there  would  be  for  /-modulation  to  relax  the  spins  anisotropically.  Unfortunately 
it  is  difficult  to  explore  this  without  much  more  efficient  techniques  for  simulating  large  spin  systems  [61, 62]. 
Extending  this  line  of  argument,  we  speculate  that  the  enhancement  arising  from  /-modulation  is  more  likely  to 
be  observed  for  a  radical  pair  in  which  the  TrpHc ,+  radical  has  been  replaced  by  a  radical  with  fewer,  smaller, 
more  isotropic  hyperfme  interactions  (e.g.  the  ascorbyl  radical  [44]).  Another  reason  in  favour  of  a  magnetically 
simpler  partner  for  the  FAD’-  radical  is  that  any  enhancement  resulting  from  /-modulation  would  be  less  likely 
to  be  obscured  by  hyperfme- induced  spin  relaxation. 

6.  Conclusions 

Previous  work  on  radical  pair  magneto  receptors  has  sought  to  identify,  and  in  some  cases  quantify,  the  factors 
that  control  their  sensitivity  to  the  direction  of  the  geomagnetic  field.  These  properties  include  the  identity  of  the 
radicals  [44, 63, 64],  their  separation  [25]  and  relative  orientation  [  5],  their  hyperfme  interactions 
[39, 44, 60, 66, 67],  and  their  chemical  kinetics  [9, 67, 68].  Clearly,  the  molecular  dynamics  of  the  radicals  must 
be  added  to  this  list  given  the  profound  effects  electron  spin  relaxation  can  have  on  the  spin- coherence.  As 
argued  here,  modulation  of  spin  interactions  by  molecular  motions  does  not  always  simply  degrade  the 
performance  of  the  sensor  by  destroying  coherence.  Under  some  circumstances,  protein  dynamics  and 
appropriate  magnetic  interactions  combine  to  produce  anisotropic  spin  relaxation  pathways  that  enhance  the 
response  of  the  compass  to  the  direction  of  the  geomagnetic  field.  Intermediate  spin  relaxation  rates  were  found 
to  be  most  effective.  Evolution,  acting  through  random  mutations  in  the  sequence  of  an  avian  cryptochrome, 
could  have  tuned  not  only  its  structural,  magnetic  and  kinetic  properties  but  also  its  dynamics  for  the  purpose  of 
magnetic  sensing.  In  this  regard,  magnetoreception  is  no  different  from  many  other  protein  functions  in  which 
dynamics  as  well  as  structure  play  an  essential  role  [69].  As  in  other  areas  of ‘quantum  biology’  [70-72], 
interactions  of  the  quantum  system  with  a  fluctuating  environment  appear  to  be  critical  for  the  quantum 
behaviour  [58, 71].  An  example,  with  a  clear  parallel  to  the  present  work,  is  ‘environment  assisted  quantum 
transport’  in  which  intermediate  dephasing  rates  enhance  the  efficiency  of  photosynthetic  energy 
transfer  [73, 74]. 
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Appendix  A.  Spin  relaxation  arising  from  /-modulation 


Equation  (7)  and  the  STD  term  in  equation  (5)  can  be  derived  as  follows.  We  start  with  the  equation  of  motion  of 
the  radical  pair  density  operator,  p  (t): 

=  —  i[H0  +  H,  it),  p(t)]  -  kp(t).  (11) 

at 

H0  is  the  time-independent  part  of  the  spin  Hamiltonian  containing  the  Zeeman  and  hyperfine  interactions  and 
the  average  exchange  interaction.  k(=  ks  =  ZcT)  is  the  rate  constant  for  the  formation  of  reaction  products  from 
the  singlet  and  triplet  states  of  the  radical  pair  (figure  1  (a)).  Hi  ( t )  contains  the  time-dependence  of  the  exchange 
interaction  arising  from  fluctuations  in  the  radical  separation: 

Hx(t)=F(t)K  (12) 

where 


Fit)  =  Hr  it))  -  /av  =  Hr  it))  -  (Jir  it))),  (13) 

hi  =  — (2Sa  •  SB  +  UJ  =  PS-  PT,  (14) 

PS  =  -E  -  SA  ■  SB  and  PT  = -E  +  SA  ■  SB.  (15) 

4  4 

/vg  zvj1  /v  yv 

P  and  P  are  the  singlet  and  triplet  projection  operators,  E  is  the  identity  operator,  and  SA  and  SB  are  the 
electron  spin  angular  momentum  operators  of  the  two  radicals.  We  assume  that  the  correlation  function  for 
/  (r  (t))  decays  exponentially  with  a  correlation  time  rc, 

(■ F(t)F(t  +  t))  =  {Fit)2)  exp  (  t/ tc)  =  ([/(r(f))  -  /av]2)  exp  (-r/rc)  (16) 

and  that  the  fluctuations  in  F  {t )  are  fast  enough  that 

(U(r(t))  -7aV]2)r2<  1.  (17) 

Then,  in  the  extreme  narrowing  limit,  a;2  pT 2  «  1(V  ft),  where  ujap  are  ff0  eigenvalue-differences, 
equation  (11)  becomes: 

+  i[Ho,  Pit)]  +  kpit)  =  -tc  {Fit)2) [hi,  [hh  pit)]]  =  -kSTDiPSpit)PT  +  PTpit)PS),  (18) 
at 

with  A:Std  given  by  equation  (7).  The  operator  part  of  equation  (18)  can  be  obtained  using  equations  (14)  and 
(19)  and 


(PS)2  =  PS; 


iPT)2  =  PT; 


PSPT  =  PTPS  =0;  Ps  +  PT  =  E. 


(19) 


We  note  that  this  approach  could  be  extended  to  somewhat  slower  fluctuations,  i.e.  outside  the  extreme 
narrowing  limit,  using  standard  Redfield  relaxation  theory. 

For  the  toy  radical  pair  discussed  in  sections  2  and  3,  the  time- independent  spin  Hamiltonian  is  given  by 
equation  (6).  More  generally,  for  a  radical  pair  with  multiple,  non-diagonal,  rhombic  hyperfine  interactions,  Hq 
has  the  form 


Ho  —  Hoa  +  Hob  —  7av^25A  •  §n  +  — hj, 

Nk 

Hok  =  7e^o  [§kx  sin  0  cos  0  +  SKy  sin  6  sin  0  +  SKz  cos  6]  +  •  A Kj  •  hp  (20) 

;= 1 

where  K  =  { A,  B } ,  and  IKj  and  A Kj  are  the  nuclear  spin  angular  momentum  operator  and  hyperfine  tensor  of 
nucleus;  in  radical  K.  NK  is  the  number  of  nuclear  spins  in  radical  K.  Note  that  the  spatial  degrees  of  freedom  are 
included  in  the  spin  Hamiltonian  by  means  of  effective  parameters  such  as  A Kj  and  /av. 

A  rough  estimate  of  ks td  may  be  obtained  by  assuming  a  Gaussian  distribution  of  separations,  centred  at 
r  =  rav,  with  standard  deviation,  a: 


Hr) 


1  1  ( 

—r=  —  eXP  I 

V27 r  a  V 


(r  ~  rav)2) 
2a2  )' 


(21) 
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Then,  using  equations  (1)  and  (7): 

J^oo 

f(r)U(r(t))  -  /av]2dr  =  4tc/o  exp ((32a2) exp (— 2/3rav)( exp  [(32a2]  -  1).  (22) 

o 

When/o  =  8  x  107 T, fl  =  14nm-1,rav  =  1.9nm,(j  =  0.1  nm,andrc  =  1  ns, equation (22) gives 
kSTD  =  2.7  x  108  s-1. 

Appendix  B.  Spin  relaxation  arising  from  D-modulation 

The  relaxation  arising  from  fluctuations  in  the  dipolar  interaction  maybe  treated  in  the  same  way  as  /- 
modulation.  The  analogous  expressions  to  equations  (13)  and  (14)  are: 

Fit)  =  D(r(t ))  -  Dav  =  D(r(f))  -  <D(r(t))),  (23) 

hi  =  -(3 SAzSBz  -  SA  •  SB)  =  -PT+  +  -PT-  -  -PTo,  (24) 

3  3  3  3 

where  P  m  is  the  projection  operator  for  the  Tm  triplet  state  (m  =  0,  db  1).  The  evolution  of  the  density  operator 
is: 


dp  it) 
d  t 


+  i[H0,  p(t)]  +  kp(t)  =  —rc  (F(t)2)  [hi,  \K  Pit)]] 


=  -  rc{[D(r(t ))  -  Dav]2)[tpSp(0PT+  +  ^Pr+p(t)PS  +  i PSp(t)Pr -  +  ^PTp(t)PS 

L  9  9  9  9 


+-psp(t)pT°  +  -pTop(t)ps  +  pr+p(t)pT°  +  pT,,p(t)PT+ 

9  9 

+pT-p(t)pT°  +  pTop(t)pT-j. 


(25) 


It  can  be  seen  from  equation  (25)  that  D-modulation  leads  not  only  to  STD  but  also  to  triplet-triplet 
dephasing;  density  matrix  elements  such  as  (T+1 1  p(t)|T0)  and  (T_i|  p(t)|T0)  (butnot  (T+1|  p(t)|T_1))are 
relaxed.  H0  in  equation  (25)  is  given  by  equations  (6)  or  (20)  with  the  exchange  term  replaced 

by  jAiv  (3Saz sBz  -  sA  •  sB). 

Using  rav  =1.9  nm,  a  =  0.1  nm,  and  rc  =  1  ns  (as  above,  equation  (21)),  the  dipolar  relaxation  rate 
constant  is: 

J^oo 

/ (r) [Dir it))  -  Dav]2dr  =  1.4  x  105  s"1.  (26) 

o 

This  value  is  substantially  smaller  than  the  corresponding  estimate  of  ks TD  for  /-modulation  (2.7  x  108  s“ :) 
because  the  dipolar  interaction  has  a  much  weaker  dependence  on  r. 


Appendix  C.  Spin  relaxation  arising  from  the  random  fields  mechanism 

A  corresponding  treatment  of  the  ‘random  fields’  relaxation  mechanism  proceeds  as  follows.  The  time- 
dependent  part  of  the  spin  Hamiltonian  is 

A(0  =  EE%w%>  (27) 


(■ FKj(t ))  =  0  and  {FKj(t)FK'j'(t  +  r))  =  (F(t)2)  exp(-r/rR), 

where  SKj  is  the;- component  (;  =  x,  y,  z)  of  the  electron  spin  operator  for  radical  KiK  =  A,  B).  Hence: 


dp  jt) 
d  t 


+  i[H0,  p(t)]  +  kp(t )  =  -tr(P(02)EE[%.  [%'.  P(t)]] 

K  j 


— — &rfr^— pit)  —  yysKj  p(t)  , 


(29) 


with 


feFR  —  2tR  {Fit)2). 


(30) 
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Appendix  D.  Spin  relaxation  arising  from  electron  hopping 

The  effects  of  the  /-modulation  caused  by  electron  hopping  between  RP0  and  Rib  (figure  1  (c))  were  modelled  by 
means  of  coupled  Liouville  equations  for  the  two  density  operators: 

— _ T0cro(t)  +  kio <J\(t)  —  /coicr0(t), 
d  t 

— 7—  =  —  kio(Ji(t)  +  /coicr0(t),  (31) 

d  t 

in  which  /c10  and  /c01  are  the  first  order  rate  constants  for  RP:  —>  RP0  and  RP0  — >  RP:  respectively  and  <j0  (0  and 
G\  ( t )  are  the  vector  representations  of  the  two  density  operators,  p0  ( t )  and  px  ( t ).  The  matrix  representations  of 
the  two  Liouvillian  superoperators  are 

L0  =  i(H0(0)  ®E-E®  H0(0)T), 

Lx  =  i  (h0(1)  0  E  -  E  0  H0(1)T)  +  y(pS0£  +  £®PS)  +  M®  E.  (32) 

The  rate  constants  k\>  and  /cf  are  defined  in  figure  1(c).  Superscript  T  denotes  transpose.  The  two  spin 
Hamiltonians  are 


Ho0)  =  -/rpo(2Sa  •  SB  + 

H0(1>  =  -/rpi  (2Sa  •  SB  +  j  +  +  H^c.  (33) 

and  H^fc  are  the  Zeeman  and  hyperfine  interactions  in  RP:  (see  equation  (6)).  It  was  assumed  that  /rp0,  the 
exchange  interaction  in  RP0,  is  so  large  that  the  Zeeman  and  hyperfine  interactions  in  RP0  have  a  negligible  effect 
on  its  spin  dynamics.  The  singlet  product  yield  was  calculated  by  assembling  a  joint  Liouville  space  for  the  two 
sites: 


so  that 


The  singlet  yield  is  then 


r*  OO 

$s  =  kbJ  Tr[Psoi(f)]df 

=  PS)eXP<-A,)^(<Q))dr 


(34) 


(35) 


(36) 


using  the  initial  condition  an  (0)  =  Ps/  Z  where  Z  is  the  total  number  of  nuclear  spin  configurations. 
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